APS/123-QED 



High Rayleigh number convection with double diffusive fingers 

E. Hage 1 and A. Tilgner 1 

Institute of Geophysics, University of Gottingen, Friedrich-Hund-Platz 1, 
37077 Gottingen, Germany 

(Dated: 17 June 2010) 

An electrodeposition cell is used to sustain a destabilizing concentration difference of 
copper ions in aqueous solution between the top and bottom boundaries of the cell. 
The resulting convecting motion is analogous to Rayleigh- Benard convection at high 
Prandtl numbers. In addition, a stabilizing temperature gradient is imposed across 
the cell. Even for thermal buoyancy two orders of magnitude smaller than chemical 
buoyancy, the presence of the weak stabilizing gradient has a profound effect on the 
convection pattern. Double diffusive fingers appear in all cases. The size of these 
fingers and the flow velocities are independent of the height of the cell, but they 
depend on the ion concentration difference between top and bottom boundaries as 
well as on the imposed temperature gradient. The scaling of the mass transport is 
compatible with previous results on double diffusive convection. 
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I. INTRODUCTION 



Double diffusive convection occurs in a ffuid whose density depends on two properties 
with different diffusion constants. The prototypical example is water whose density depends 
on temperature and salt concentration. This particular combination is relevant in oceanog- 
raphy and the interest in double diffusive convection historically originated in this field. 
More applications have since emerged in geophysics, astrophysics, and engineering^. Double 
diffusive convection differs from ordinary Rayleigh-Benard convection in that two fundamen- 
tally new flow structures can appear: Layers and fingers. Fingers are vertically oriented long 
and narrow regions of up- or downwelling fluid which can occur if the fluid property with 
the large diffusion coefficient (for example temperature) imposes a stabilizing gradient on 
the fluid, whereas the property with the small diffusion constant (the salt concentration in 
the oceanographic example) is unstably stratified. The basic principles of finger formation 
are well understood since linear stability analysis already predicts convection with fingers 2 . 
Layers dividing the fluid volume into several sections in the vertical direction with little 
advective transport between the layers is frequently observed in double diffusive systems, 
irrespective of which component is stabilizing and which is destabilizing. The necessary 
conditions for layer formation remain controversial^-. 

The experiments reported here deal with a double diffusive convective system which 
forms fingers but no layers so that various scaling laws obeyed in the finger regime can be 
determined accurately. 

A popular experimental setup for the study of double diffusive convection has been water 
with gradients of salt and sugar concentrations. The diffusion coefficients of salt and sugar 
differ by a factor of about 3. However, the ratio of the diffusivities of temperature and salt 
is much larger, typically around 100. The large aspect ratio of fingers makes their numerical 
simulation difficult. Experiments on the other hand have a problem in maintaining a steady 
state. Many experiments start from an initial distribution of salt and temperature and let 
the salt stratification disappear in the course of time. Some results may then depend on 
whether an experiment was started from a continuous stratification or a step function in 
the salt distribution-. One is therefore interested in experiments capable of keeping both 
constant temperature and concentration differences between top and bottom boundaries. 
Maintaining a temperature difference is simple, but the salt stratification is problematic. One 
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solution is to use an apparatus with permeable membranes as top and bottom boundaries^. 
Water tanks with constant temperature and salt concentration placed outside the convecting 
volume on the other side of each membrane guarantee steady conditions in the experimental 
cell. 

An alternative solution to the problem consists in using an electrochemical system. Con- 
vection in electrochemical cells with vertical electrodes has been studied at least since 
1949 9,10 . Even an application to double diffusion with horizontal gradients appeared^. 
Electrochemical convection was later noted in a cell with horizontal electrodes^ and devel- 
oped into an analogy for Rayleigh-Benard convection^. The principle of the system used 
most frequently and also in this paper is to fill a cell with a solution of CuSO^ and to 
place copper electrodes at the top and bottom of the cell. When an electric current is sent 
through the cell, copper dissolves from one electrode and is deposited on the other electrode. 
A non-uniform spatial distribution of copper ion concentration is responsible for buoyancy 
and convection. In order for copper ions to behave in the same way as salt or temperature 
in a Rayleigh-Benard experiment, one needs to make sure that the ions diffuse and are ad- 
vected, but do not experience a force in the electrical field due to the potential difference 
applied between the two electrodes. This is achieved by dissolving a large concentration of 
another electrolyte such as H 2 SO^ which does not participate in the chemical reaction at the 
electrodes and which screens the electrical field in the bulk of the cell. The electrical field 
separates ions of the electrolyte which accumulate in charged layers of microscopic thickness 
next to each electrode so that the copper ions move essentially in a field free environment. 
Refs. k] and 14 provide more details about the electrochemical system and its relation to 
Rayleigh-Benard convection. 

In the experiments presented here, the two copper electrodes are regulated in temper- 
ature, with a cold cathode at the bottom and a warm anode at the top. Temperature is 
thus stabilizing and the ion concentration is destabilizing. Mass transport is conveniently 
measured in this system because it is directly related to the current flowing through the 
electrodes. Since the working fluid is transparent, we also use particle image velocimetry 
(PIV) to measure velocity and finger sizes. We are able to vary the chemical and thermal 
Rayleigh numbers over several orders of magnitude. It will be shown below that the thermal 
stratification is important for the mass transport and leads to fingers even if the thermal 
Rayleigh number is much smaller than the chemical one. The main results will be finger 



widths, flow velocities and mass transport in convection with double diffusive fingers in a 
statistically steady state. 

The next section presents the experimental apparatus and methods. The results are 
summarized in table |T] and discussed in the third section. 



II. EXPERIMENTAL PROCEDURES 

This section will describe a novel experimental realization of a double diffusive convec- 
tion cell in which temperature and ion concentration are the two diffusers. It is useful to 
first introduce some nomenclature and definitions to describe the system. Double diffusive 
systems are characterized by four control parameters: There are two diffusivity ratios, the 
Prandtl number Pr and the Schmidt number Sc, 

Pr= V - , Sc=^ (1) 

in which v stands for the kinematic viscosity of the convecting fluid, k for its thermal 
diffusivity, and D for the ion diffusion constant. There are furthermore the thermal and 
chemical Rayleigh numbers, Rclt and Ra c , 

gaATL 3 gflAcL? 

Ra T = , Ra c = — (2) 

KV Dv 

with g the gravitational acceleration and L the cell height. The two expansion coefficients a 
and (3 determine variations of density p around a reference state with density, temperature, 
concentration and pressure po, T , Co and po via 

ru \ / co,po,po ru \ / To,po,po 

Both a and are positive. The sign convention for the Rayleigh numbers used in this paper 
is such that negative Rayleigh numbers indicate a stable stratification. This implies that 

AT = Tbottom ~ T t op , Ac = C top — Cbottom (4) 

with the subscripts indicating the boundary at which temperature T or concentration c are 
evaluated. Another dimensionless number in common use which will become useful below 
is the density ratio A: 

Rcit k a AT 
= ~R^D = JA^' ( ' 



A quantifies the ratio of thermal and chemical buoyancy forces. 
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The electrochemical aspects of the experiment are essentially the same as in ref. 
cell was made of two copper electrodes of 1cm thickness painted with varnish on all surfaces 
which were not in contact with the electrolyte. Water from thermostats was circulated 
through pipes welded onto the copper plates in order to regulate their temperature. The 
sidewalls were made of four plexiglass plates 1cm thick epoxied together to form a rectangular 
frame. Four of these frames have been used, all with a cross section of 20cm x 20cm and 
heights of L = 2cm, 4cm, 8cm and 20cm. 

The cell was filled with a solution of CuSO^ in 1 molar sulfuric acid. The solution was 
prepared by dissolving CuSO^ ■ {H 2 0)^ at a few tens of mmoljl. In between measurements, 
samples of the solution were extracted from the cell in order to check with optical absorption 
spectroscopy that the copper concentration remained constant. For PIV measurements, 



polyamide particles 100/zm in diameter were added to the solution. The materia 
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constants 



necessary for the determination of Pr, Sc, Rclt and Ra c were taken from ref. 

The concentration difference Ac is not directly under control in these experiments. It is 
known only for suitable settings of the electrolytic cell. If no current flows through the cell, 
Ac = and the concentration is everywhere equal to the average concentration cq. At small 
voltages applied to the electrodes, the current rises for increasing voltage. Only diffusion 
transports the copper ions close to the boundaries so that a concentration gradient must exist 
near each electrode. The copper ions are depleted near the cathode and their concentration 
is increased near the anode. The so called limiting current is reached when the copper ion 
concentration at the cathode is down to zero^. It is concluded from symmetry that the 
anode concentration is then 2 • cq, and therefore the concentration difference across the cell 
is Ac = 2 • Co- Increasing the voltage further cannot increase the current any more and a 
plateau is reached in the characteristic curve of current versus voltage. An increase beyond 
the limiting current is only possible for voltages large enough to allow another chemical 
reaction, such as water dissociation. 

The voltage range for which the limiting current is obtained depends on the cell height, the 
concentration cq, and AT. Preliminary measurements of the voltage characteristic for AT = 
—4 K and AT = were performed for all heights and concentrations in order to determine 
the appropriate voltage used later to create and maintain the concentration difference of 
Ac = 2 ■ Cq. 
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FIG. 1. Sh as a function of Ra c in a fluid with uniform temperature. The dashed line indicates 
Sh = 0.14S'c~ ' 03 i?a^' 297 , which is the functional dependence proposed in id. The continuous line 
is given by Sh = 0.052Ra° c u . 

The Sherwood number, which is the same thing as the chemical Nusselt number, is 
directly proportional to the number of ions transported from top to bottom divided by the 
purely diffusive current, so that the Sherwood number can be determined by 

Sh = jL A (6) 
zFDAc K ! 

if j is the current density, z the valence of the ion (z = 2 for Cu 2+ ) and F Faraday's con- 
stant. A transference number, which must sometimes be taken into account when computing 
Sherwood numbers^ is so small in the present case that it was neglected. 

The double diffusion experiment was preceeded by another experiment in which the 
CuSO^ solution was kept isothermal in cubes of side length 5cm, 10cm and 15cm. The 
dependence of Sh on Ra c obtained from this preliminary experiment is shown in fig. [TJ 
This figure contains also one point at Ra c w 10 9 obtained in one of the cells later used for 
the double diffusive experiments (the one of height 2cm). These results should be identical 
to measurements of the Nusselt number in a thermal convection experiment. Figure [1] also 
shows the Sh(Ra c ) dependence extrapolated from measurements in ref. [lg, performed in 
a cylindrical cell at Rayleigh numbers between 2 x 10 7 and 3 x 10 10 and Prandtl numbers 
between 4.3 and 1350. In the interval of Rayleigh numbers where both experiments overlap, 
there is a discrepancy of about 20%. This may be due to the different cell geometries, but 
discrepancies of this order of magnitude are not unusual when comparing different Rayleigh- 



FIG. 2. Shadowgraph picture taken across the cell of height 80mm 15min after the voltage was 
switched on for Ra c = 1.15 x 10 11 and Rclt = —9.70 x 10 7 . The bright horizontal bar is an artefact 
of the illumination. The picture covers the entire height of the cell and a width of 85mm. 

Benard experiments at these Rayleigh numbers^. 

A typical experimental run of double diffusion convection was performed as follows: Care 
was taken when filling the container that no air was trapped inside the cell. Then the top 
and bottom plates were set to the desired temperatures and enough time was allowed for a 
linear temperature gradient to establish. The voltage was then applied instantaneously and 
the current observed. During transients lasting from a few minutes to one hour, the cur- 
rent oscillated with decreasing amplitude and period until the limiting current was reached. 
After that time, the current fluctuated by less than 1% and a statistically stationary finger 
pattern filled the whole cell. This pattern could be observed in shadowgraph pictures taken 
in parallel, an example of which is shown in figure |2j The precise nature of the initial con- 
ditions turned out to be irrelevant for the final state. A few runs were started from uniform 
temperature and copper concentration, and the water circulation through the thermostats 
and the voltage were switched on simultaneously. This procedure led to longer transients 
but did not change the measurements thereafter. 

PIV measurements were performed in the stationary state. The light beam of a 50mJ 
Nd:YAG laser system with two IR laser heads illuminated a vertical plane with a width of 



FIG. 3. Velocity field obtained from PIV measurements at Ra c = 1.05 x 10 9 and Rut = —3.81 x 10 5 
(left panel) and Ra c = 1.14 x 10 9 without any stabilizing temperature gradient (right panel). The 
picture shows the entire cell height of 20mm and extends over a width of 12mm (left) or 45mm 
(right). The lines (red online) indicate direction and velocity of the flow, the shade of gray depends 
on the vertical component of velocity which varies from — 0.2mm/ s to 0.2mm/s (left) and from 
— 0.1mm/ s to 0.7mm/ s (right) in going from dark to bright. The left panel shows the average of 
60 pictures taken during an interval of 15s in order to reduce the noise contained in any single PIV 
velocity field, and the right panel shows an average over 30 pictures. 

approximately 3mm in the middle of the cell. The time separation between the two pulses 
varied between 50 and 1000ms depending on the Rayleigh numbers. A camera took pictures 
of an approximately square region extending from the top to the bottom electrode. An 
example of a velocity field deduced from PIV measurements is shown in figure [3J Once 
the fingers were established, they remained surprisingly steady, keeping their positions for 
30mm or longer. For each velocity field the rms velocity was calculated and averaged over 
all pictures to give the velocity V. From there, the non-dimensional measure of the velocity, 
the Reynolds number Re, is computed as 

Re = Y±. (7) 

For the determination of the finger size, the vertical velocities in the PIV velocity fields were 
Fourier transformed along the horizontal and the spectra averaged over all horizontal lines 
and all pictures. The location of the maximum of the spectrum yields the finger thickness. 
Independently, the finger thickness could simply be measured on the pictures with a ruler. 

The accuracy of the determination of d is limited by the number of fingers visible on any 
one picture. The relative error on d is typically 10%. The error on V depends on the setting 
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of the PIV system and could be as large as 20 %. The limiting current is the most accurate 
of our measurements with an error of a few percent. There is also an uncertainty on the 
control parameters. The measurements at small AT suffer most from fluctuations in the 
temperature regulation. Moreover, the material constants are based on extrapolations from 
relatively few measurements^ which adds up to an uncertainty on the Rayleigh numbers on 
the order of 10 %. 
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TABLE I. Summary of experimental results. The first three columns give the cell height, the 
applied temperature gradient, and the concentration of CuSO^. The next four columns con- 
tain the nondimensional control parameters (Prandtl, Schmidt, thermal Rayleigh and chemical 
Rayleigh numbers) and the last three columns the nondimensional measured quantities (finger 
width, Reynolds and Sherwood numbers). 



III. RESULTS AND DISCUSSION 



The main quantitative results are contained in table HI The most important qualitative 
observation from shadowgraphs and PIV is that fingers exist at all. This is remarkable 
because apart from two exceptions, |A| < 1 in the experiments and the total density is 
unstably stratified so that fingers are not necessary for convection to start. One might also 
think that at |A| = 10" 2 , the presence of the weak stable thermal stratification should be 
irrelevant for the chemical convection, but it is not. Fingers still appear. They are of course 
absent for A = 0. The isothermal experiments in figured] have A = and it could be verified 
at the occasion of those preliminary experiments that a single convection roll forms in a 
cubic cell in this case and that the aspect ratio of convection rolls in the cell of height 20mm 
is compatible with what is known from ordinary Rayleigh-Benard convection^. 

Some more information about the structure of the fingers can be deduced form table [H 
The finger thickness always exceeds the concentration boundary layer thickness A, which 
can be computed from the measured Sherwood number as 

A is small compared with the cell height L in all cases. Since the fingers extend across the 
entire height of the cell, they carry fluid directly from one boundary layer to the other. In a 
horizontal cross section, they must therefore consist of a core of size A, which is the detached 
boundary layer, surrounded by the entraining fluid filling a cross section of size d. The ion 
exchange between neighboring fingers is small as long as the distance over which ions diffuse 
during the time it takes to transit from one boundary to the other is small compared with 
the finger size, i.e. as long as (DL/V) 1 ^ 2 jd = (Re Sc)~ l l 2 L/ d « 1. It can be seen from 
figure H] that chemical exchange between fingers can be neglected for the experiments listed 
in table HI and that this will not be true any more around |A| w 5. 

If there is no significant chemical diffusion across fingers, it is a simple matter to estimate 
Sh for Sh » 1 from the advective transport: 

Sh^Sh-l = ^-Vc (9) 
DAc v ' 

where c denotes the concentration anomaly in a finger averaged over its cross section, c = 
\(\/d) 2 Ac in a finger with circular or square cross section. This assumption is not compatible 
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A AT/ L = -10.00 K/cm 

4 AT / L = -5.00 K/cm 

T AT / L = -3.00 K/cm 

► AT / L = -2.00 K/cm 

• AT / L = -1 .00 K/cm 

♦ AT / L = -0.75 K/cm 
■ AT / L = -0.50 K/cm 

★ AT / L = -0.25 K/cm 

* AT/ L = -0.10 K/cm 
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FIG. 4. (.Re Sc)- l / 2 L/d as a function of |A|. The symbol shapes indicate the applied temperature 
gradient, the colors the height of the cell: 20mm (red), 40mm (dark blue), 80mm (green) and 
200mm (light blue). 



with the data in table [I] If on the other hand fingers form sheets, c 
shows that 

1 



Sh = -(Re Sc L/d) 



1/2 



\\/dAc. Figure [5] 



(10) 



is a good representation of the data from which we conclude that fingers have a lamellar 
shape. 

Pr and Sc vary by less than 30% in table |T] so that a dependence on these parameters 
cannot be extracted from the data. It must be kept in mind that the prefactors in the power 
laws given below potentially depend on both Pr and Sc. 

Let us now turn to the scaling of the finger width. Which parameters do we expect 
to determine this width? Long and narrow fingers should not be affected by far away 
boundaries, d should thus be independent of L. The fluid inside fingers carries the ion 
concentration from the boundary layer the finger started from. There are little diffusive losses 
during the transit, so that d should depend on Ac. Temperature on the other hand diffuses 
much more rapidly and fingers somewhere in the bulk do not know about the boundary 
temperatures. Temperature must enter the expression for d through the vertical gradient, 
AT/L. A thickness d which depends on Ac and AT/L but is independent of L translates in 
nondimensional terms into a law of the form d/L oc {Raxl 11 Raj 2 with 471 + 372 = — 1. The 
scaling known from linear stability analysis, d/L oc |i?ay|~ 1//4 belongs to this family of power 
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(Sc Re L/d) 1 ' 2 

FIG. 5. Sh as a function of Re Sc L/d. The line is given by eq. (fTU|) . The symbols are the same 
as in figure HI 




FIG. 6. d/L as a function of |i?a-r| 1 ^ 3 i?ac . The line is given by eq. (jllj) . The symbols are the 
same as in figure HI 

laws. However, the measured finger sizes in table U depend systematically on Ra c . 71 and 
72 have been determined from a linear regression applied to \og(d/L) = A + 71 log l-Ra^l + 
72 log i?a c . Our best fit obtained with that procedure yields 71 = —0.32 and 72 = 0.086, 
which is close to 

y = O.^lRarl-^Ray 9 . (11) 

The quality of this fit can be judged from figure [6j The important point is that the exponents 
obey 471 + 372 = — 1 so that eq. (TXTT) is compatible with the general picture of fingers 
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FIG. 7. Re as a function of |i?ay| l / 2 Ra c . The line is given by eq. (I12p . The symbols are the 
same as in figured 

delineated above. 

The velocity should in the same fashion depend on Ac, AT/L, but not on L. This means 
for the nondimensional variables that Re oc I-Rot] 73 Raj 4 with 473 + 374 = 1. An additional 
argument can be made concerning the dependence on the gravitational acceleration g. If 
buoyancy is balanced by viscous friction, one has to find v oc g. However, a balance of linear 
terms is questionable because the linear terms all occur in linear stability analyses predicting 
fingers growing without bounds as a function of time. Saturation of finger velocity at finite 
amplitudes must involve the nonlinear terms, even if they are small. A balance between 
buoyancy and the advection term leads to v oc g 1 ^ 2 . A well known example of such a 
situation is Rayleigh-Benard convection with a Rayleigh number near its critical value -Ra cr ; t . 
The velocity near onset is small so that the Reynolds number is small, but experiments and 
weakly nonlinear analysis 19 show that velocity is proportional to (Ra — Ra^it) 1 ^ 2 oc g 1 ^ 2 . 
Returning to double diffusion, the additional requirement v oc g 1 ^ 2 fully determines the 
exponents 73 and 74 to be 73 = —1/2 and 74 = 1. And indeed, a good fit to the data of 
table fl] is close to (see figure Ej) 

Re = 10- 6 \Ra T \- 1/2 Ra c (12) 

A linear regression applied to the logarithm of eq. (|T2l) yields as best fit Re = 2.9 x 
10~ 6 \Rarp\~°' 3S Ra^ S7 . However, Re is our measurement with the largest scatter and eq. (fl2"|) 
is only a marginally worse fit than the result of the linear regression. Another reason to 
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FIG. 8. Sh as a function of \Rcit\ 1 ^ 12 Rat . The line is given by eq. (|T3l) . The symbols are the 
same as in figured 

prefer eq. (TT2I) over the direct fit is that eqs. (TTUl [TT| [T2|) now fix the dependence of on 
i?a^ and Ra c to be 

S/i = 0.016|i?a T |- 1 / 12 J Ra;y 9 (13) 

which is shown in figure |8] and which is also recovered from a direct fit by linear regression 
to the logarithms, whose result is Sh = 0.017|i?aT|~°' 095 -Ra°' 45 . Since Sh oc L/\ and since 
we may again argue that A ought to depend on Ac, AT/L but not on L, we expect a scaling 
Sh oc \RaT\ ls Ral 6 with 475 + 376 = 1. The exponents in eq. ( fTBl obey this constraint. 

According to eq. ( IT3l) . the Sherwood number for finite \Rclt\ should become larger than 
the Sherwood number for zero temperature stratification shown in figure [U if Ra c exceeds 
\Rclt\ by five orders of magnitude. Such a combination of parameters has not been reached in 
the experiments. It is not a priori impossible that the Sherwood number for finger convection 
becomes larger than the Sherwood number for convection without fingers, but it seems more 
likely that the corresponding combinations of control parameters mark the limit of validity 
of eq. ( TI~3"1) . Such a limit must exist somewhere because Sh — > 00 for \Rclt\ — > in eq. ( |T3l) . 

We finally compare the results of this section with previous work. The scalings of Sh, Re 
and d/L show that fingers in our experiments take the form of sheets. On the theoretical 
side, as analytic treatments of the planform selection problem became increasingly realistic, 
the predictions went from square cells to rolls and back to square cells^^i. Experimentally, 
both square cells^. and cells in the form of rolls have already been observed^. 
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Some experiments have determined a variation of finger widths with \Rar\ compatible 
with \Rclt\~ 1 ^ — This is the size of the fastest growing instability in a fluid layer with 
uniform temperature and salt gradients. Starting from a step change in temperature and 

-1 / 3 , the same exponent 
22l The scaling for the flow 



salinity, ref. |6| found numerically a finger width scaling as \Rar 



that appears in eq. ( II ip and which is also compatible with ref. 
velocity, eq. (|T2|) has apparently never been observed before. 

Past measurements of mass transport in finger convection^— have found good agreement 
between experimental data and a relation of the form Sh oc Ra}J z f(Pr, Sc, A) with an 
undetermined function /, which is sometimes called the 4/3-law because this relation states 
that the mass flux varies as Ac 4//3 at constant Pr, Sc and A. Arguments in favor of this type 
of relation come from dimensional^ and asymptotic analysis^. Eq. (fl3|) can be rewritten as 
Sh oc Ra 1 J' i \A\^ 1 ^\Rax\ 1 ^ m ■ The variation in |i?ay| 1//36 is much too weak to be detectable 
in the experiment, so that eq. ({TBI must be considered to be indistinguishable from the 
4/3-law. 



IV. CONCLUSION 

With an electrochemical technique, double diffusive convection can be investigated with 
a tabletop experiment for Pr f» 9, Sc ~ 2200 and chemical Rayleigh numbers ranging 
from 5 x 10 8 to 10 12 while varying the density ratio from 10~ 2 to 1. Fingers appear in 
all circumstances. The main quantities related to mass transport, the finger thickness, the 
flow velocity, and the concentration boundary layer thickness, are all independent of the 
cell height L, and depend on the concentration difference between top and bottom Ac as 
well as on the imposed temperature gradient AT/L, but not on the temperature difference 
AT alone. The scaling for the Sherwood number proposed in eq. (fTBl is experimentally 
indistinguishable from the classical 4/3-law. 

Perhaps the most remarkable finding is that fingers occur in the first place, considering 
the small stabilizing temperature gradient used here. It remains an open question whether 
there is an abrupt or a continuous transition from fingers to ordinary convection rolls which 
are at least as wide as they are tall. Eq. fTTTj) remains correct down to the weakest stable 
stratification controllable in the experiment. In the case of a continuous transition, eq. (TTT]) 
predicts that the small thermal Rayleigh number of \R(It\ ~ Ra l J 3 is necessary for d w L. 
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This indicates that conditions suitable for finger formation are much more widespread than 
linear stability analysis suggests 2 -. 
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